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The spectral line shapes for hydrogen-like heavy ion emitters embedded in strongly correlated 
two-component electron-ion plasmas are investigated with numerical simulations. For that purpose 
the microfield fluctuations are calculated by molecular dynamics simulations where short range 
quantum effects are taken into account by using a regularized Coulomb potential for the electron- 
ion interaction. The microfield fluctuations are used as input in a numerical solution of the time- 
dependent Schrodinger equation for the radiating electron. In distinction to the standard impact and 
quasistatic approximations the method presented here allows to account for the correlations between 
plasma ions and electrons. The shapes of the Ly a -line in Al are investigated in the intermediate 
regime. The calculations are in good agreement with experiments on the Ly a -line in laser generated 
plasmas. 
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I. INTRODUCTION 

Measurements of emission and absorption spectra of atoms and ions are one of the most important tools in plasma 
diagnostics [H-0 ■ They allow for investigations of the properties of various laboratory and space plasmas. In particular, 
spectral line shapes may be analyzed to yield a wealth of information on the plasma parameters provided, however, 
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(-h ■ that the data are compared to accurate computations of the spectral line broadening. 

Traditionally, in the theory of spectral line broadening in plasmas ion effects in most cases were calculated within 
the quasistatic approximation, while the electron perturbation was believed to satisfy the impact approximation. The 
description of this method (Standard Theory (ST)) is given, e.g., in Refs. [H-Q- The separation of perturbations into 
ion and electron parts, in general, cannot be made without a loss of accuracy, although it is argued Q that in many 
^T) ' cases it is justified. A more serious problem, however, is that each of the ion and electron parts often needs to be 
considered beyond the limits of the quasistatic and impact approximations. Inparticular, the ion motion in plasma 
leads to the so-called ion dynamics effects. It was first shown theoretically 0, @ and soon found in experiments 
that the ion dynamics can be responsible for significant corrections to the spectral line widths. In order to 
advance the calculations beyond the ST, several numerical methods have been developed. Among the first is the 
method developed in Ref. [lOj], where a computer code was used to simulate the ion motion along straight paths, 
while the electron contribution was calculated using the impact approximation. The method was further improved in 
Ref. [ll[ by using molecular-dynamics (MP) simulations for the ions, thus accounting for interactions between the 
, radiators and the ion perturbers. In Refs. (l2| - [l5| the motion of both ions and electrons was numerically simulated. 
The particle motion was simulated using straight path trajectories, which is applicable when the correlations between 
the perturbers and radiator are neglected. Later, the area of applicability was extended by using h ype rbolic [l6j and, 
recently, exact paths for the perturbers (see, e.g., Refs. [l7l - [T9j and references therein). In Refs. (171 [l8j the model 
was based on one-component plasmas (OCPs) treating the full spectrum as a superposition of the electronic and 
ionic contributions and thus neglects the influence of the attractive interactions between electrons and ions. This is 
well justified for weakly coupled plasma where the ionic and electronic fields can be handled separately for spectra 
modeling. But there is an increasing number of experiments of interest which are far beyond such parameter regimes 
(see, e.g., Refs. (20l - [23j ). In such cases a simple superposition of the electronic and ionic fields becomes insufficient due 
to strong nonlinear effects and the total field in a two-component plasmas (TCP) should serve as the starting point 
for spectra modeling. The related microfield distribution (MFD) and the line shapes including the full attractive 
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electron-ion interaction has thus been attracting more and more attention and has already been studied, e.g., in 

Refs. mm. 

The present paper is the continuation of Ref. [l8[ but we treat here the ions and electrons on an equal footing by 
concentrating on the TCP. For that purpose we perform MD simulations which span the entire range between the 
impact and the quasistatic approximations. We solve the time-dependent Schrodinger equation for the radiator in the 
fluctuating microheld generated by the plasma particles. The MD simulations in conjunction with the Schrodinger 
equation to study the line shapes in a TCP have been previously considered in Ref. (19J. Our model is thus similar 
to that considered in .19] but, in addition, it allows to treat the interaction of the emitted photons with radiating 
electron. 



II. MD SIMULATIONS 



A two-component classical plasma of electrons and ions (with the charge Ze) is in equilibrium with the temperature 
T completely described by the coupling parameters r a( g with a,j3 = e,i. Introducing the mean electron-electron (a e ), 
ion-ion (a^) and electron-ion (a) distances through the relations, a~ 3 = 47TOq,/3, a -3 = 4wn/3 (where n = n e + m is 
the total plasma density with n e = Zrii) these parameters are defined as 
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Here q e = —1, qi = Z, e% = e 2 /A-keq and £o is the permittivity of the vacuum. 

It is well known [32|, l33j that, to avoid the collapse of the classical system of electrons and ions, the Coulomb 
electron-ion interaction potential must be replaced by the pseudopotential with a regularized short-range behavior. 
In this paper we consider the electron-ion pair interaction potential — e|g^u e i(r), where /3 = i,R, que is the charge 
of the radiator (throughout this paper the index R refers to the radiators) and 

u ei (r) = ^(l-e- r /% (2) 

which is regularized at small distances. The cutoff parameter 5 may be qualitatively thought of as a classical emulation 
of the electron thermal de Broglie length. For large distances r > 5 the potential becomes Coulomb, while for r < 5 the 
Coulomb sin gularity is removed and u ei (0) — 1/5. By this the short range effects based on the uncertainty principle 
are included |3ll - [33l |. 

For a classical description of a plasma the electron degeneracy parameter e , i.e., the ratio of the thermal energy 
and the Fermi energy must fulfill O e = ksT/Ep > 1. Since an ion is much heavier than an electron this condition 
is usually fulfilled for ions. Therefore one can expect that the regularization given by Eq. (j2J is less important for 
ions than for electrons. Furthermore, scattering of any two particles is classical for impact parameters that are large 
compared to the de Broglie wavelengths % a p = {2irh 2 / ' fia/sksT) 1 / 2 , where /j, a p is the reduced mass of the particles a 
and j5. Typical impact parameters are given by the Landau lengths, Alo.0 — e 2 s \q a qp\/kBT . Its ratio to the de Broglie 
wavelengths is given by 

^ei\q a qp\ a 

O a {3 = 7} T ■ (3) 

Zj A a /3 

Note that a ee < a e i <C uu- A classical description of the scattering events in the TCP is valid if a ee > 1. 

A collective length scale is given by the Debye screening radius, for a TCP Ad = a/ (Sr^) 1 / 2 . The plasma frequencies 
for electrons and ions uj pa = (iTrq^naeg /rria) 1 ^ 2 with a = e,i set the collective time scales uj~^ and u;^ 1 for electronic 
and ionic subspecies, respectively. Due to their large mass ratio the electrons and the ions move on very different 
timescales. Moreover, for a nonrelativistic treatment the thermal energy of the particles must be smaller than their 
rest energy. Since this is important only for electrons we require that ksT <C m e c 2 . Also the validity of the dipole 
approximation for plasma-radiator interaction used in Sec. IIIII requires that the characteristic length scale of the 
plasma microfield must be larger than the effective atomic length scale az = as/Z, where as is the Bohr radius. 
Since this length is ~ a the dipole approximation is valid when a > az which is usually fulfilled for heavy radiators. In 
this paper we consider hydrogen-like ions as radiators in a completely ionized TCP, i.e. we assume that qn = Zr — 1, 
where Zne is the charge of the nucleus of the radiating ion. 

Let us now briefly discuss the limitations of the MD model arising from the classical description of the electrons, 
i.e., from the neglect of the quantum effects in the short-range electron-ion interaction. As shown in Ref. (3~ij two 
constraints on the parameter 5 determining the regularized potential ([2]) must be considered. In the parameter regime 
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FIG. 1: Normalized electric microfield distributions for Al 13+ TCPs with k B T = 500 eV, n e = n , F ee = 0.037, T i4 = 2.65 (left) 
and n e = 2no, P ee = 0.046, Tu = 3.34 (right) as a function of the electric field in units of Eh, Eq. (UJ, for different values of 6. 
The lines with and without symbols represent the MFD from the MD simulations and PMFEX model, respectively. The open 
circles are the MFD obtained from the folding of an electronic and an ionic OCP, see Eq. ©. The Holtsmark distributions for 
a TCP (see Ref. [2q|l are shown as dot-dashed and dotted lines. (Online color: www.cpp-journal.org). 



when a significant fraction of the simulated electrons is found in the quasibound states the simulated Z and Zr are 
effectively reduced and the MD simulations are not adequate. Thus, the parameter S must be chosen large enough 
to suppress the formation of the classical bound states of electrons. On the other hand the condition S < a must 
be fulfilled so not to affect the free electron density at r ~ a. The probability to found an electron within a volume 
r < S from an ion is estimated by W ~ (4ir5 3 /3)n e g , e ,(0), where the electron-ion radial distribution functi on Q ei(r) 
can be approximated by the nonlinear Debye-Huckel expression g e i(r) ~ exp [f3 e Z e 2 s u e i{r)] (see, e.g., Refs. [251 |3l|) 
with /3 e = l/k^T. Note that in the case when Z » 1, the capture of an electron on the quasibound orbit reduce 
the effective Z and Y e i for the next one, so a significant fraction of electrons can stay free even for T e i<i/5 3> 1. The 
minimal value of W(S) occurs at S — T e ia/3. It is thus clear that for T e i > 0.5 the capture of the electrons onto 
classical orbits becomes important and the significant contribution from the classically bound electrons cannot be 
avoided in the MD models with point-like classical electrons. 

The electric microfield distribution (MFD) P(E) plays a central role for the line shape. Models for this distribution 
exist in the limits of an ideal plasma |35| , a weakly coupled plasma (3|| and for very strongly coupled plasmas [13] ■ 
For intermediate cases an effective independe nt -p article model known as Adjustable Parameter Exponential (APEX) 
approximation has been developed in Refs. [38l44fj| | for an ionic OCP. It rests essentially on the pair distribution 
function and has been tested by comparison with MD and Monte-Carlo simulations. Recently in Refs. [25Tj27j we have 
suggested the theoretical models named PMFEX (Potential of Mean Force Exponential approximation) and PMFEX+ 
which turn out to be a very reliable approaches for calculating the MFD of a TCP with attractive interaction. In 
order to cover the entire range from small to large plasma parameters we use here classical MD simulations which 
have been described in detail in Ref. [4l| (see also Refs. (25l - [27| ). As an example the normalized MFDs from PMFEX 
and MD are compared in Fig. 1 where the electric microfields are scaled in units of the Holtsmark field Eh for a TCP 
(see [25[ for details) 
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with an effective charge Z and C = (87T/25) 1 / 3 . These distributions were obtained from ensembles of fields taken at 
a charged reference point which is chosen to be one of the plasma ions. 

The MFDs for Al 13+ TCP with temperature 500 eV and with coupling parameters T ee = 0.037, T u = 2.65 and 
T ee = 0.046, Tu — 3.34 are shown in left and right panels of Fig. Q] respectively, for different values of S. The 
density of plasma electrons is measured in units of no = 5 x 10 23 cm~ 3 . The dot-dashed and dotted curves are, for 
comparison, the Holtsmark MFDs for a TCP with regularized Coulomb potential. Here these MFDs depend on Z 
and 5 as discussed in [25[. To demonstrate the importance of the attractive interactions we also plotted the MFDs 
Pq{E) resulting from the corresponding electronic and ionic OCPs with T ee and F^, respectively (open circles). To 
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that end the distribution Qo(E) = Pq(E) / (An E 2 ) of the total field E = Ei + E 2 is calculated as 

Qo(E) = J <m 1( m 2 5 (E — Ei — E 2 ) Q e (E 1 )Q i (E 2 ) (5) 

from the MFD of the ionic OCP at a charged point Qi(E 2 ) and of the electronic OCP at a neutral reference point 
Q e (Ei). The distribution Qq(E) thus represents the MFD in a TCP assuming that the ion-electron attractive 
interaction is absent. Here Q e (Ei) and Qi(E 2 ) are taken from MD simulations of an OCP. As the thermal motion 
of the particles is suppressed with increasing coupling the distributions P(E) and the mean electric fields are shifted 
towards smaller values as shown in Fig. [TJ This figure also shows the importance of the attractive interactions in 
plasmas. The behavior of the MFD with respect to the variation of the parameter 8 is particularly noteworthy. For 
fixed coupling parameters the maximum of the MFD shifts only slightly to lower field strengths E with increasing 
5 (Fig. [T}, while the maximum itself increases with 8. This is related to the largest possible single-particle field 
|i? e (0)| = e/(8irso8 2 ), which an electron can produce at the ion. Thus the nearest neighbor electronic MFD vanishes 
for electric fields larger than |i? e (0)|, and smaller 8 will result in larger contributions to P(E) at higher fields E with 
a corresponding reduction of P(E) at small fields. Therefore, the formation of the tails in the MFD and enhancement 
of the electric microfield at small 8 may have important influence on the spectral line shapes of the radiating particles. 
Further examples for charged and neutral radiators , togeth er with a detailed discussion of the limits of the PMFEX 
treatment at increasing coupling, are given in Refs. [25l427| . 



III. WAVE EQUATION FOR A RADIATOR 



In this section we describe the solution of the wave equation for a hydrogen-like ion coupled to the time-dependent 
electric microfield. The microfield fluctuations in the plasma are calculated by the MD simulations as discussed in 
SecHU 

We consider a hydrogen-like ion in a time-dependent electric microfield. The Hamiltonian is the sum of Ho 
describing the unperturbed ion and a dipole term i?j n t = er • E(i) for the interaction between the bound electron 
(distance r from nucleus) and the microfield E(£), 

H = H + Ant ■ (6) 

The electron moves in the potential of a nucleus with charge Zrc In the present application it turns out that it 
suffices to start from the non-relativistic Schrodinger equation 



for the time-independent electronic state \a) with energy E a — huj a . Here p is the momentum operator and a is 
a multiindex including radial, angular momentum and spin quantum numbers. The present calculations are done in 
the configuration space corresponding to the solutions of Eq. ([7]) . In order to discretize the continuum a boundary 
condition (r|a) = is imposed at a radius r = Rq, which is chosen sufficiently large in order to avoid an influence on 
the final results. The radial wave functions with this boundary condition are still confluent hypergeometric functions, 
but the radial quantum numbers of bound states are not integers any more [13] • In order to obtain a finite basis 
the (former) continuum states are cut off at sufficiently large quantum numbers. Alternatively the continuum could 
be handled by forming wave packets with a width that must be adjusted appropriately (42[. In Refs. [13, EH the 
time-dependent equation H^(t) = ih^(t) with the Hamiltonian © has also been solved on a grid for the electron 
wave function [H, H3] . This is more advantageous for the description of the continuum and it is easier to implement 
the interactions between the radiator and the plasma particles beyond the dipole term in Eq. ([6]) . On the other hand 
spatially extended states require very large simulation boxes. Quite generally in the present context the solution on 
the grid is more expensive numerically than working in configuration space. We adopted the latter for the subsequent 
calculations. 

At high Zr relativistic corrections must be considered and also the spin should be treated as a dynamical variable. 
It turns out that in the cases considered here it suffices to include the first-order fine-structure shift [45| 

ZWs\E n \ ( 1 3 



ABn , = -r«™ , (8) 

n \7y An) 

where 7y = j + 1/2 and 7y = 1 at I ^ 1 and I = 0, respectively. Here n is the principal quantum number of the 
hydrogen-like ion, E n = —Z^EsIn 2 is the corresponding energy (Eb = e 2 s /2aB is the Bohr energy), j = Z ± 1/2 is 
the total angular momentum quantum number and as ~ 1/137 is the Sommerfeld constant. 
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The dipole interaction er • E(t) between the radiator and the plasma is time-dependent and possibly strong. Going 
beyond the second order treatment ol Rcf. 4(| we use the interaction picture with the unperturbed basis states 
given by Eq. ([7])- The time-dependent Schrodinger equation with the total Hamiltonian (j6|) can be solved using 
Dirac's method. The perturbed electron wave function is represented as a sum of wave functions of the unperturbed 
Hamiltonian with time-dependent coefficients c a (t) 

*(t)=53^(*)e- <u - t |a> ■ (9) 

OL 

A substitution of Eq. @ into the time-dependent Schrodinger equation and orthogonality of the spatial wave functions, 
i.e., (a = 5 a p, gives the set of coupled ordinary and linear differential equations 

t a (t) = -Iffi (t) • £ e^cp® (a |r| /3) (10) 



which is solved iteratively. Here hjj a p, is the transition energy between atomic states a and j3, i.e., Lo a p = ui a — uip. 
Within the dipole approximation the transition rate per unit time and energy interval I(ui) for the emission of photons 
is proportional to the power spectrum of the dipole operator (47j . It is defined as the square of the absolute value of 
the Fourier transform of the expectation value of the dipole operator. Hence we introduce 



1(c) = 



2e|w 3 ' 



with r — > oo and the expectation value of the dipole moment 



f &{t)e iut dt 
Jo 



(11) 



d(i) = £e^*c^)c a (iH/3|r|a). (12) 

Let us now consider a transition a — > g downwards to a state \g) which is nearly filled, i.e. cp(t) = 5p g . Then the 
dipole moment in Eq. (|12l) is calculated with respect to the state \g) and 



V(ff|r|a) f c a (t)e^-^dt 
Jo 



2 

(13) 



Our numerical model includes also the interaction H C1 of the radiating electron with the emitted photons which, 
however, has been neglected in Eq. ([5]). In this case as shown in Refs. [l7], [H| the total radiated power given in 
Eq. (|13|) is underestimated for the excited radiators where \c g \ 2 < 1. This can be compensated by dividing through 
the time-averaged occupation probability of the lower state 

lM-> JM 2 ■ (14) 
(\c g (t)\)t 

The subsequent calculations will be done in this dipole power spectrum approximation (DPSA). We have tested the 
validity of this approximation in the wide range of plasma and radiator parameters by comparing explicitly l(oS) with 
the spectrum obtained from the total Hamiltonian H' = Hq + H ln t + H ei . We have found that in the parameter 
regime considered in Sec. IIVI the DPSA is justified as the emission of radiation through the interaction H C1 changes 
the occupation probabilities of the radiator's states on a much slower scale than the fluctuating electric microficlds. 
However, our preliminary results show that the electron-photon interaction H ei may have an important contribution 
to the wings of the line especially in the case of light emitters. We intend to take up further studies on this issue in 
a separate paper. 

At this stage we have neglected the feedback of the radiator's excitation to the plasma. In this respect the plasma 
particles move as if they had an infinite mass. As they have a finite velocity it appears as if the radiating electron 
is embedded in a plasma of infinite temperature. Accordingly the time evolution of the total system will lead to an 
equal population of all electronic states. As the time-dependent feedback could be implemented only at a very great 
expense in the MD simulations we enforce a canonical equilibrium state of the plasma and the radiating electron by 
modifying the interaction in Eq. (|10[) according to 

er • E(f) -> e- p ' 6o ^er ■ E(t)e^«*°/ 2 . (15) 

The time-dependent equation (JTDJ) describing the coupling of the microfield to the radiator is then solved for an 
ensemble of typically thirty independent microfields which yields the mean emission as well as the statistical error. 
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TABLE I: Doppler broadening (FWHM) of the Ly a -line in an Al plasma. 



k B T (eV) 10 10 2 5 x 10 2 10 3 10 4 10 5 
FWHM (eV) 0.08 0.26 0.57 0.81 2.57 8.12 



IV. RESULTS 



Using the theoretical background introduced so far we present in this Section calculations of the shape of the Ly a - 
line of Al 12+ radiating ion embedded in a A1 13+ -TCP in a wide range of plasma parameters. As we have mentioned 
in Sec. Mil for heavy ions the relativistic corrections, i.e. the fine structure of the levels must be accounted for. Using 
the fine structure shift in Eq. (J5J the unperturbed Ly Q transition energy becomes 



ftuj 



L, a = l Z 2 R E B (l + C j Z 2 R a 2 s ) 



with Ci/2 = |g, C 3 / 2 = Jq anQ Ior Ar ~ ' lons nLd Ly a — if^o-i ev ana rawLy Q 
respectively. 

We start from the line as it is broadened by the Al 13+ ions and electrons in the plasma. Then we fold with the 
weighted fine structure shift Eq. §8$ which is /iAw ~ 1.29 eV according to 



jg and for Al 12+ ions ^>Ly Q 



1728.1 eV and fiw 



(16) 

1729.4 eV with j = 1/2 and j = 3/2, 
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(17) 



and account for the Doppler effect [T], Q which broadens a line with the unperturbed frequency uj according to a 
Gaussian distribution 



2? (Aw) 



1 



2tt<j 



exp 



where Aw = w — ujq and 
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(18) 



(19) 



Here fkvw is the full width at half maximum (FWHM) of the line and Mr is the radiating ion mass. Some values of 
FWHM huw for aluminum are shown in Table [I] Finally, we fold to account for the experimental resolution. 

We discuss now the simulated Ly Q -line profiles at solid state densities no ^ n e ^ 4no and at k B T = 500 eV. Some 
results for the Ly^-line shape without Doppler broadening and LS coupling are shown in Figs. [5][SJ In all cases the 
ratio az/a <C 1 is small and the use of the dipole approximation in Eq. ([5]) is fulfilled. The line broadening and shift 
towards lower photon energies (redshift) is clearly visible in Fig. [5J which shows the line profile at fixed temperature 
500 eV and different densities. Here the regularization parameter S = A e j is determined as the thermal wavelength. 
With increasing density the influence of the plasma effects on the line shape becomes more pronounced. At very 
large densities from n e — 2.6n (dashed line) and up to the value n e — 4n (dash-dotted line) there is hardly any 
broadening and the line is only redshifted towards lower photon energies by the plasma effects. In this high density 
regime the fluctuating electric fields become sufficiently strong to cause asymmetric shapes due to nonlinear coupling. 

To gain more insight we now fix the plasma temperature (500 eV) and the density (n e = n ) and show in Fig. [3] the 
line profile for different regularization parameters S (the lines without symbols), 5 — 0.08a (solid line), 6 = 0.1a (dashed 
line) and S = 0.4a (dotted line). Here the thermal wavelengths of the electrons are chosen as the relevant lengths 8 at 
which a smoothing of the ion-electron interaction due to quantum diffraction becomes effective. For comparison we 
also calculate the line shapes for non-isothermic plasma with different electronic and ionic temperatures T e = 500 eV 
and Ti = 50 eV, respectively (the lines with symbols). As shown in Fig. [3] the width of the lines decreases with 
increasing parameter <5, i.e. by "softening" of the ion-electron interaction. Besides, keeping the electron temperature 
unchanged and decreasing the ionic temperature leads to an additional broadening of the lines and this is visible for 
more Coulomb-like interactions with small 5. 

In Fig. 2] we put together the results obtained for two values of plasma densities (no and 2no) and the parameter 
d. In the top and bottom panels we take S = l.baz and S = l.Qaz, respectively, where az is the effective Bohr 
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FIG. 2: Simulated Ly a -spectra of a Al 12+ radiating ion embedded in a A1 13+ -TCP of a temperature of 500 eV and solid state 
densities n e = no (solid line), n e = 1.6no (dotted line), n e = 2.6no (dashed line) and n e = 4no (dash-dotted line). The spectra 
are normalized to the area under the curves. (Online color:www. cpp-journal.org). 




FIG. 3: Simulated Ly a -spectra of a Al 12+ radiating ion embedded in a 

A1 i3+_ TCp of a density of 

n e =no. The lines without 

and with symbols correspond to the equilibrium (T e — Ti — 500 eV) and non-equilibrium (T e — 500 eV, Ti = 50 eV) TCPs, 
respectively. The regularization parameter is 8 = 0.08a (solid line), 8 — 0.1a (dashed lines) and 8 = 0.4a (dotted lines). The 
spectra are normalized to the area under the curves. (Online color:www. cpp-journal.org). 



radius of Al 12+ . Note that for two plasma densities no and 2ng the chosen values of 6 in units of az are equivalent 
to 0.08a, 0.1a and 0.4a, 0.5a, respectively, in units of the Wigner-Seitz radius a of a TCP. We also demonstrate the 
influence of the electron-ion attractive interaction on the spectral line shapes plotting the spectra Tq(w) resulting from 
a superposition of the electronic and ionic OCPs (dashed lines). Io(u>) is calculated by folding the spectra I & {uS) and 
Xi(oj) which are obtained from simulations of the radiative transitions of a radiator embedded in an electronic OCP 
and of the ionic OCP, respectively. The microfields E(t) in an electronic and ionic OCPs are simulated at a neutral 
and charged reference points, respectively. The spectrum Xo(ui) thus represents the line shape in a TCP assuming 
that the ion-electron attractive interaction is switched off. As shown in Fig. [4] the width of the line now turns out 
to be highly sensitive to both the choice of 8 and the density n e , with a much stronger dependence on n e for smaller 
5. The influence of the high-electric field tails in the MFDs at small 5 on the spectral lines is now clearly shown in 
Fig. ID Smaller 5 results in higher electric fields which broaden the spectral lines and reduce the peak intensity. More 
precisely we observed that the line width behave approximately as HAuj ~ n e /S. 

In the following we will compare our calculations with the results of experiments performed in Garching [22[ where 
an Al plasma is created by the irradiation of the target with laser pulses of 150 fs duration at an intensity of a 
few 10 17 W/cm 2 . The systematic investigations carried out in Refs. [U [H| snow that the standard (quasistatic 
and impact) approximations become doubtful if the plasma density reaches that of the solid state. In the last years 
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FIG. 4: Simulated Ly a -spectra of a Al 12+ radiating ion embedded in a A1 13+ -TCP (solid lines) of a temperature of 500 eV 
and solid state densities n e = no (left panels) and n e = 2no (right panels). The spectra are normalized to the area under 
the curves. The parameter 5 is 1.5az and 7.6az in the top and bottom panels, respectively. The dashed lines are the spectra 
obtained from the folding of an electronic and an ionic OCP. (Online colonwww.cpp-journal.org) . 

experiments have approached this regime, see, e.g., (2(j| - [22l |. We note that the theoretical model discussed so far 
assume a homogeneous equilibrium plasma. Obviously this is not the state in which the laser leaves the target 
after the irradiating pulse. In particular self-absorption due to plasma inhomogeneities leads to an additional line 
broadening which is difficult to analyze. Fortunately there has been considerable experimental progress to reduce the 
self-absorption [22j. 

Earlier experiments on the Ly Q -line in Al 12+ at solid state density [2(| [2l[ were subject to self-absorption in the 
cooler and less dense surface regions of the target. This can be prevented by using thin (to reduce absorption) target 
layers with sharp boundaries (to enhance homogeneity). For that purpose a 25 nm Al target layer was embedded 
in solid carbon at depths ranging from d = 25 nm to d — 400 nm [22| . With increasing depth the expansion of the 
Al layer is suppressed and the homogeneity of the Al plasma is improved. In Fig. [5] we compare our simulations 
with the experimental results (filled circles) for d — 400 nm and k^T — 500 eV from which the underground has 
been subtracted. The fine structure and the Doppler broadening are taken into account as described above. Then the 
simulated Ly Q -lines are folded with the experimental resolution (0.9 eV, FWHM) and compared with the experimental 
line assuming densities 5 x 10 23 cm -3 (dashed line) and 10 24 cm -3 (thin solid line). At these two densities the plasma 
parameters for the A1 13+ -TCP are Tu = 2.65, r ee = 0.04 and T zi = 3.34, r ee = 0.05, respectively. All curves in Fig. [5] 
are normalized to the peak intensity. Finally, the position of the simulated line must be redshifted by 2 eV. This is 
the dense plasma line shift (DPLS) Ref. [l[ due to the screening of the electron-nucleus interaction by hot background 
electrons. Assuming a Debye-screened interaction instead of the r _1 -Coulomb potential first-order perturbation 
theory yields a shift of the required magnitude. A comparison of the two simulated curves in Fig. [5] allows to conclude 
that the remaining uncertainty in the determination of the density of the target is of order 10 23 cm~ 3 . Our results 
show that the quantum mechanics of close electron-ion collisions is important over and above the plasma redshift. If 
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FIG. 5: Comparison of the experimental line [52] (filled circles) with our simulation results (dashed and solid lines), i.e. with the 
solid curves from the top panels of Fig. [4] after taking into account the Doppler broadening, the LS-coupling, the experimental 
resolution and a redshift (see text). The experimental line is subtracted by the underground. Here the curves are normalized 
to the peak intensity. (Online color:www. cpp-journal.org). 

the quantum diffraction parameter S is fixed at physically reasonable values near the effective Bohr radius az, our 
calculations favor a somewhat larger density than no as proposed in Ref. [22}. Clearly a more quantum mechanical 
treatment of the electron component in the plasmas is desirable. 

V. CONCLUSIONS 

In this paper we have presented a model for spectral lines that works without some assumptions which underly the 
conventional impact and quasi-static approximations. In particular we (i) consider two-component plasma (TCP) 
with attractive interactions between electrons and ions, (ii) account for the strong Coulomb correlations between 
plasma particles, (iii) account for radiator states including the continuum, which are not directly involved in the 
transition, (iv) allow for a non-perturbative treatment. 

We have compared our model with recent experiments on Al targets and found good agreement for the Ly Q - 
transitions. The more exact treatments beyond the standard approximations will become highly desirable in connec- 
tion with experiments at higher densities and temperatures at the planned (X)FEL facilities. 

A critical discussion of our results suggests further improvements, (i) The dipole approximation in Eq. © for 
the interaction of the microfield with the radiator suffices for the present experiments |22j. In even denser plasmas 
one must account for close collisions between the radiator and the plasma particles with a quadrupole term in the 
expansion of the interaction and finally with an exact treatment [46j |. (ii) Relativistic and spin effects beyond the 
simple fine structure given by Eq. ([8]) can be taken into account by treating the radiator with the Dirac equation 
p2l ]. (iii) The major He-like satellite is well separated from the Ly Q -line in the experiment (22[. However, there 
will be closer satellites due to spectator electrons in higher configurations, which may affect the "red" shoulder of 
the line. For spectators in the continuum this effect merges into the DPLS. The satellites impose a challenge as 
they offer an additional tool to determine the temperature of the plasma, see, e.g. Ref. 0. For that purpose one 
has to solve the multi-electron wave equation, for example in the relativistic case the Dirac equation [48|. (iv) At 
high densities of the laser-produced plasmas the self-absorption is likely to occur. In general, the consideration of 
radiative transfer in dense plasmas is involved. To estimate the effect of the self-absorption in the calculation of the 
line profile one-dimensional monolayer model can be used which has already been successfully applied in Ref. (49[ 
(see also references therein). 
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